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The acoustic Faraday rotation in the 4f paramagnet Tb3Ga50i2 has recently been observed by 
Sytcheva et al (arXiv: 1006 .0141 ). As in earlier examples the rotation angle per unit length ol 
transverse acoustic modes was found to depend linearly on sound frequency. Existing theories 
for this effect consistently require that it should vary with the square of the frequency. In the 
present work a solution for this long-standing problem is provided. We propose a model based on 
magnetoelastic interactions with 4f quadrupole moments that includes both acoustic and optical 
phonons. The symmetry allows a direct and induced coupling between the latter. This leads to an 
indirect acoustic Faraday rotation via the field induced splitting of doubly degenerate optical phonon 
modes. It varies linearly with frequency in accordance with experiment and dominates the rotation 
angle. It also explains the observed resonance of the rotation angle in the field range between 17-20 
T. The mechanism is of general validity for non-Bravais lattices and applies to previous examples 
of the acoustic Faraday effect. 

PACS numbers; 72.55.-|-s, 73.50.Rb, 71.70.Ch 

I. INTRODUCTION 

In crystals with a fourfold uniaxial or higher symmetry transverse acoustic phonons which propagate along axis 
direction are twofold degenerate with perpendicular polarizations of the displacement vector u. Equivalently they 
may be described by left (L)- and right (R) handed circularly polarized phonons which belong to complex conjugate 
representations of the group Gq of the phonon wave vector q which is aligned with an axis. The degeneracy of these 
complex conjugate modes is ensured by time reversal invariance. If the latter is broken by application of a magnetic 
field the degeneracy may be lifted. Depending on field direction this leads to the Faraday rotation for H || q or 
Cotton-Mouton effect for H ± q in the long-wave length limit. Only the former will be considered here. This means 
that for fixed sound wave frequency oja the wave numbers ql and qr of circular polarised modes will be different. In 
an ultrasonic experiment a linearly polarised transverse mode is generated which may be described as a superposition 
of circular modes with equal amplitudes and phase shift tt/2. Since the latter have different wave numbers this means 
that the polarisation direction of the linearly polarized mode will be continuously rotated as it propagates along the 
crystal axis. The displacement field of the acoustic wave is then given by 

u(z, t) = A{'x. cos (f) + y sin cf)) exp(za;at — iqz) (1) 

Where A is the amplitude, z the propagation (axis) and x,y the polarization directions. Furthermore average wave 
number q and Faraday rotation angle 4> are given by 

9=^(9L+gfl) (j) ^ ^{qL - qR.)z = (j)' z (2) 

Here 4>'{H) is the Faraday rotation angle per unit length. 

The splitting between q^ and qn wave numbers due to time reversal symmetry breaking will only be effective 
when the sound waves can couple strongly to magnetic degrees of freedom. Therefore the acoustic Faraday effect was 
first observed in magnetically ordered 3d compounds such as yttrium-iron garnet (YIG) [l|, [2[ and in paramagnets 
with magnetic impurities [3|-|5|- In the former it is due to magnetoelastic coupling caused by the strain dependence 
of the anisotropy energy of ordered moments Q. It has also been proposed for S=l/2 paramagnets and for 
paramagnetic 4f compounds with general CEF split level scheme [sj. In the latter case the Faraday effect is due 
to the magnetoelastic coupling to (2J+1) crystalline electric field (CEF) split 4f states where J is the total angular 
momentum. This mechanism was first found in the paramagnetic phase of CeAl2 [9|. In a magnetic field H | 
q the magnetoelastic coupling leads to nondiagonal quadrupolar susceptibilities of the 4f CEF states resulting in 
a Faraday rotation of the sound wave polarisation. Both mechanisms lead to a Faraday rotation angle (j)' which 
increases quadratically with frequency uja- However experimentally in CeAl2 the rotation angle turned out to increase 
linearly with sound frequency. This discrepancy has been unresolved sofar. Recent Faraday rotation experiments on 
paramagnetic Tb3Ga50i2 or terbium-gallium garnet (TGG) (Ref. ITol ) gave the same linear frequency dependence. 
Therefore it is clear that the existing theories need to be extended to explain these experimental findings. We note that 
the Faraday rotation angle in the magnetically ordered systems like YIG was also predicted to vary with the square 
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of the sound frequency [l|, |3| . However this has never been tested experimentally because only a single frequency has 
been used in existing experiments. 

II. MAGNETOELASTIC COUPLING MECHANISM FOR ACOUSTIC AND OPTICAL PHONONS 

The displacement field of sound waves may be described by elastic strains which couple to the CEF quadrupolar 
(or higher order multipolar) moments of the 4f shell. This causes a temperature dependent rcnormalisation of sound 
velocities (elastic constants) [ll| and under suitable conditions a Faraday rotation of the polarisation 8] due to the 
splitting of complex conjugate (R,L) acoustic modes. It is natural to expect a similar effect for doubly degenerate 
complex conjugate optical phonons. This has indeed been found in Ref. [l3 with Raman scattering in a magnetic 
field and explained in Ref. Il3l . An interesting possibility arises in non-Bravais lattices with an atomic basis: In the 
case that acoustic and optical phonons for q along an axis cause distortions of the same symmetry they will be 
coupled, both directly and indirectly via quadrupole excitations of 4f CEF states. Such an a-o coupling will lead 
to an additional contribution to the Faraday rotation which has not been considered before. This new contribution 
turns out to solve the problem of the frequency dependence of the rotation angle, both for TGG and the previous 
case of CeAl2. 

Our starting Hamiltonian therefore has to include a term (Hph) describing both acoustic and optical phonons. In 
addition it contains a part (i?4/) corresponding to the CEF- split 4f electrons and a phonon-4f electron interaction 
part Hph^if. The coupled system is then described by 

H = Hph + + Hph-4f (3) 

First we turn to the phonon part. Restricting to modes with q along z they may be written in terms of acoustic 
(•^q/i = '^q^ + '^-q/i) ^^'^ Optical ($q^ = ^qp + A^q^) phonou coordluatcs where jj, — x,y denotes the two polarisations 
and CLq^, are the usual phonon creation operators. 

= E ^^q'^q^e''''''; Q,i^) = E Qq^qM^'"^- (4) 
q q 

Here Ri denotes the lattice sites, furthermore we defined (N,M = number and mass of unit cell, respectively). 

Q-^^iqi2MNoj-^)-i; Q° = (2M7Vc^°)-^ (5) 
The total phonon Hamiltonian is then given by 

Hph = E + ^] + II <[4p^qA< + ^] + II ^q*qM<M (6) 

qM q^ q^ 

Here we defined Aq — iq{E/4M){uj^u}^)^^/'^ where E is the acoustic-optical (a-o) coupling constant between long 
wavelength acoustic strains and internal displacements of the optical mode. Note that in this representation Hph 
has not been completely diagonalised since there is an interaction term reflecting the fact that acoustic and optical 
phonons belong to the same type of representations of the group Gq of the wave vector. 

The lattice vibration modes of the garnet structure RE3AI5O12 which is isostructural to TGG (space group 0\^) 
have been classified and discussed in Ref. [13. Due to the large unit cell there are 240 vibrational modes for q ^> 0. 
In addition to the acoustic Ti„ mode (one longitudinal, two transverse) there are 17 Ti„ and 14 Tig optical modes 
and in addition 14 T2g and 16 optical modes. For finite q || [001] the symmetry is reduced from cubic to uniaxial 
. Therefore the triply degenerate Ti and T2 modes split into singlet A and doublet E modes for both acoustic and 
optical phonons. The acoustic E modes correspond to the sound waves with (x,y) polarisation for propagation along 
z. Because of their E symmetry they will couple with all doubly degenerate optical E modes through the background 
force constants and the true acoustic eigenmodes are complicated superpositions of long wavelength and internal 
displacements. Therefore in a model based on magnetoelastic interactions it is preferable to start with modes having 
pure acoustic strain and internal displacements and include the a-o coupling explicitly as in Eq. Only in this 
case the magnetoelastic interaction can be written in terms of the simple strain and displacement amplitudes given 
in Eq. ([5]). In the model we will include only one representative optical mode. This may be taken as the lowest Tiu 
or T2g mode which is around 100 cm~^ in RE3AI5O12 [li) and presumably also in TGG. 
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The 4f part of the Hamiltonian including the effect of an apphed magnetic field is given by 

H^f = ^ Ea\T an) {^anli " h ^ Jzi^) — ^ ^^an |f q„ ) (f |i (7) 
ani i ani 

Where Ea and \Tan) are the zero-field CEF level energies and states respectively and h = gjiiBH{J = 6,gj = |) is 
the field variable. The last expression in the above equation is written in terms of Zeeman split CEF energies Ean 
and states iFcm). 

The lattice vibrations couple to the 4f states by changing the CEF potential due to induced local distortions [Tlj . 
In this way the E- type displacements of acoustic and optical modes will couple to E-type quadrupolar moments 
iOxz,Oyz) of the 4f shell. Then the phonon-4f or magnetoelastic Hamiltonian is given by 

where fi = x,y denotes polarisation for phonons and the quadrupolar operators (//=xz,yz,) are defined by 

Oa^z = JxJz + JzJx = ^ [{J+Jz + JzJ+) + {J-Jz + JzJ-)] 
Oyz = JyJz + JzJy = Yi ^^•^+'^' + •^'•^+'^ ~ ("^-^^ + '^"'^-)] 

with J± — Jx iiJy. Furthermore ga and go are the magnetoelastic coupling constants for acoustic (a) and optical (o) 
modes respectively. 



III. DYSON EQUATIONS FOR THE PHONON PROPAGATORS 

To obtain the propagating modes in the presence of magnetic ions and external field we have to set up Dyson's 
equations for the phonon propagators of doubly degenerate a,o modes. As a first step, for better understanding of 
the a-o coupling term in Eq. ^ we treat the case without 4f-phonon coupling and magnetic field. Then the phonon 
propagator is simply given by 

where the diagonal elements are the unperturbed inverse phonon propagators for a,o modes. The nondiagonal terms 
mix the bare modes to the true vibrational modes whose frequencies are obtained from det{D~^) = or 

This leads to renormalised dispersions for both modes. Since we may assume ^ w", at least in the long wave 
length limit the above equation leads to 

|2 ,,,o |x 12 



^ = ^ni_2M!<i; ^° = ^°ri + 2M<i (12) 

q qL ,,o2 ,,,aJ' q qi ,,o2 ,,,oJ ^ ' 



-q -q ^q 

This shows the repulsion of a and o modes due to their interaction. In the long wavelength limit where uj'^ = Vaq and 
Wq = Wo — Dq^ we obtain a renormalised sound velocity for the acoustic mode 

Va=Va[l-'^]=Va[l-\-&] (13) 
Wo ^ 

where 7^ = {E / AM)"^ / {vl^iOo) and E^ = 47^/a;o is the dimensionless a-o coupling constant. The complementary 

,0 
q 



correction to the optical phonon dispersion may be expressed as w° Wo — Dcf' where 



,14, 
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Thus the a-o couphng of acoustic strains and internal strains characterised by the constant E leads to a reduction of 
sound velocity and a reduced dispersion of the optical mode. 

Now we treat the full problem with magnetoelastic interactions included. We will show that the a-o coupling (E) 
plays an essential role in the acoustic Faraday rotation mechanism. For finite field and ga,5o 7^ the two polarisations 
have to be included explicitly for both modes. Then the total phonon propagator may be written as a 4 x 4 matrix 

i^-(q,-)-(^Xt ^^?) (15) 

where the diagonal 2x2 blocks are given by the following expressions (the momentum label q has been omitted for 
simplicity) 



where s = a, o is the mode index. Furthermore the non-diagonal block describing a-o coupling is given by 

AM^f/l -Sh{u^)\ (17) 



with A = *(|A| + Sd)- This matrix is anti-Hermitian with = —A. The diagonal {S^) and nondiagonal (S*^) magne- 
toelastic self energies for a given mode (s = a,o) and Sh, Sd for the a-o coupling are obtained as (sl [l3|: 



Sdii^) = ^K^o)*5'((0..0..))L; Shiio) = -{cuaLUoy^fiiO^^Oy,))'^ (18) 

Here we defined g — (gago)^ and Oq,^ = Oq^ — {Oap)- The double brackets denote the dynamical susceptibility which 
are explicitly given in Sect. IIVI Furthermore we used the definitions 

N\9aQl? = \vaq~gl or ~gl- 



2 ^" CaV, 

Arl /00|2 ^ ~2 ~2 Ho 

N\goQ^\ = T^t^ogo or .g„ - 



(19) 

where Ca = pv^ is the elastic constant (in the cubic garnets = C44) with p — M /Vc being the mass density (M, Vc— 
unit cell mass and volume, respectively). 

The new phonon frequencies modified by magnetoelastic coupling as well as the a-o coupling are then obtained 
by finding solutions of defD~^(q, w) = 0. For that purpose it is convenient to transform the matrix propagator to 
circular polarised phonon coordinates L(-l-) and R(-) according to = — i^y) and 0_ = + i4>y) for 

acoustic and similar for <^± circular optical modes. Then we can regroup the 4x4 propagator in pairs of {(p-,^-) 
and (0+, $+) modes which leads to the decoupled 2x2 circular mode propagators 



1 (^2_~±2) _A± 



where the renormalised circular mode frequencies are given by (s=a,o) 

= w^q ± 2a;^qS'^(w) 



and the effective a-o coupling is obtained as 



-sq 

'^'q = - -^Stiiu)] (21) 



A± - i[\X^\+Sd{uj)±S,M] (22) 
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The last two terms in this equation are the dynamical correction caused by magnetoelastic interactions to the static 
background a-o coupling Aq. The total remaining a-o coupling in Eq. (pHl) may be easily diagonalised by finding the 
zeroes of the determinant of the inverse 2x2 propagators. This leads to the secular equations 

_ ^±2)(^2 _ ~±2) ^ 4^^^^^^|A±|2 (23) 

For the acoustic modes we may approximate uia — ^aq ^ <^oq- This leads to the renormalised circular mode frequencies 

= ^Iq ± 2a.a,5;jKq) - 4a.,q^oq|A±|2-l^ (24) 

The equivalent relations for the circular optical modes will not be considered further here. In the following we neglect 
the last term in Eq. ([2^ because for a; Sh{i^) ^ oj and therefore it is negligible compared to Sd{L^)- 



IV. FARADAY ROTATION ANGLE 



From the above result we may easily compute the Faraday rotation angle. In the Faraday configuration the frequency 
LOa is fixed and L,R modes have different wavenumbers qL,qR which may be derived from Eq. (|24p as 



I 
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(25) 



where the renormalised sound velocities of L,R acoustic modes are obtained in the limit ujaq — > as 



Va 
Et 



_ -2/'-, S^^^a 

= =^;a[l-5a((0..0..))^,] 



(26) 



and the split optical mode frequencies are given by Eq. (f2T|) . Here we defined the dimensionless bare a-o coupling 
constant by i? = {EaEo)^ with Ea = E/{caVc) and Eo — E/{Muj1) which is equivalent to E"^ = 4j'^/lUo (see 
below Eq. (llSC l.We note again that Et contains the dynamical magnetoelastic corrections to the bare a-o coupling 
{g — (gago)^)- These corrections may be strongly field dependent as discussed in Sect.|Vl Under the approximation 
(f)' — ^{ql ~ Qr) ^ ^(^i + 9r) — 9 we obtain the acoustic Faraday rotation 



HI" 



— 9 

Va 



(27) 



Using the explicit form of the quantities appearing in this expression and 9 = ^ = x '^here A denotes the ultrasonic 
wave length we get the final result 



A \Va 



1 



~9i{{O.^Oy.)):^+Ef-{-:^-—^ 



(28) 



with split optical phonon frequencies and renormalised sound velocities explicitly given by 



■,±2 



- l + ~gl{{O.A.))l^±~gl{{O^.Oy,))'^ 



1 



(29) 



The Faraday rotation in Eq. ((28)) consists of two parts. The first one , g^) which was already discussed in 

Ref. d is due to the direct coupling of long-wavelength strains to the quadrupole moments of the 4f shell. The second 
(j)'^ (~ Ef) and indirect part derived in the present work is due to the coupling of long-wavelength strains to the 
internal optical displacements which in turn couple to the 4f quadrupoles. The latter leads to a splitting of circular 
polarised optical phonons and subsequently (via Et) to a splitting of circular acoustic modes. Therefore an indirect 
optical phonon contribution 0^ to the acoustic Faraday rotation appears. We note that the common prefactor j 
in Eq. (|28l) is linear in the sound frequency. Since the first term in the bracket is also ~ uia (see Eq. ([50]) ) one has 
4''a ^'^^ direct acoustic contribution. This is because for g — > the acoustic mode frequencies and therefore 



6 



also the splitting of +,- (L,R) acoustic modes has to vanish. This contributes the additional frequency factor in 4>'^. 
This is different for the indirect part Because the sphtting of optical ujf modes stays finite for q ^ the second 
term in the brackets of Eq. psp is simply a (field dependent) constant and one finally gets (j)'^ ^ Ua- For small Wq 
the second part will always dominate and one has 4>' ^ LOa for the total rotation angle. This solves the long-standing 
puzzle of the frequency dependence of acoustic Faraday rotation in paramagnetic 4f compounds. 

In the following we give more explicit forms for the Faraday rotation by using the expression for the quadrupolar 
susceptibilities [8, 13] 

((0«0,))L = E' ^'"^^' x^^^^ ^i" + {kT)-'Y.Pi{im^){l\dp\l)5^, 

Im ^ ^Im I 

{{OAs)): - Y^'^^^f^A^ (30) 

Here Kim — Ei—Em are the excitation energies between CEF levels in an external field and pi = e^'^'^'''^ / e^'^"/*'"^ 
are their occupation numbers. The prime denotes summation over terms with Ei ^ Em- Furthermore Oq = Oa — {Oa) 
and the (anti-) symmetrised matrix elements are defined by 

S'rp = \[{imm){m\O0\l) + {m\0^\l){l\Op\m)] 

A% - \mOM){MO^l) - {rn\0^m\Of,\Ta)\ (31) 

It is implied here that |^), |to) denote the eigenstates in the external field. In our case a,ji — xz^ yz and the notation 
will be simplified to A^zyz = A^y and Sxzxz = Sxx etc. Expanding the inverse of the split optical frequencies 
in the magnetoelastic coupling cf^ in Eq. (|28p and inserting the expressions for quadrupolar susceptibilities we obtain 
after some algebra the Faraday rotation angle (f)' = + (j)'^ as 



t(— \9a}^ — ; T^Axy + [—] [Etgo) > — ; x^^^y 



Remembering that j uja this formula shows explicitly that the first part ip'^ ~ because we may take — 
in the denominator: A typical sound frequency of I'a = '^a/2vr = 0.5GHz corresponds to = 0.15 K. This energy 
is much smaller than all other energy scales involved, in particular smaller than the lowest CEF excitation energy 
At = 57.3 K (see sect.|V| and even smaller than the thermal energy at the measuring temperature 1.4 K which gives 
the scale of the CEF broadening. Therefore the above approximation is justified and the iOa frequency dependence 
of the first part is simply given by {iT/X)u!a ^ ^a- the second part due to indirect coupling to optical phonons 
the acoustic frequency appears only through the prefactor (tt/A) leading to a linear behaviour (p'^ ^ uja ■ One may 
estimate under which condition the linear (indirect optical) part dominates the quadratic (direct acoustic) part of 

the rotation angle. This happens when g^uja (s^) iEtgo)^i^o- Using Wo — "^o and assuming ga — go,Et ~ E one 

obtains the condition E ^ {uja/u>o)^ ■ For the typical uJa/^-n = 0.5 GHz and lOq — 150 K this means we must have 
E :s> 0.9-10"^. Assuming from Eq. ([T^ that one has a modest change of sound velocity {va — Va)/va — 0.5-10"^ due to 
the bare a-o coupling only we get E = 0.1, which is much bigger than the required value. Therefore we conclude that 
under any noticable a-o coupling the indirect (~ uia) contribution to the Faraday rotation will dominate the direct 
acoustic one (~ w^) and the sum ,i.e., the total observed Faraday rotation will be linear in uta to a good approximation. 

The fact that optical phonon frequency Wo and sound velocity Va are renormalised by the coupling to the diagonal 
quadrupolar susceptibilities enters only in a non-essential way through the modification of prefactors in Eq. (|32p . 
Nevertheless we give their explicit expressions for completeness. The renormalised optical phonon frequency is obtained 
from the selfconsistent solution of 



~ 2 2 



'AimiPl -Pm) ^Ir, 



'^Tr^r' s^ (33) 

Im 



and the renormalised sound velocity from 

E't - 3^(E'%^^- + (^^)"' T.p^{i\o.Ai){i\dxz\i)) ~ {E~go)'Y^ ^':f' 'J--^ 
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FIG. 1: Zeeman splitting of simplified cubic r3(0) — r5(At) — r2(A3 = 3.3At) level system. Only transitions between Fa — F5 
states corresponding to full lines have nonzero quadrupole matrix elements (Eqs. (|A3IA4|) ') and contributions to the quadrupole 
susceptibility. The empirical matrix elements used for the Zeeman splitting are m^^ = —4.5, 77135 = —2.0, mig = 0.0. Field 
scale: H = 28.6 T for h/At=0.5. The Fj - FJ excitation ( A^") crosses the lowest optical phonon mode {ujo ~ 2. 5 At) around 
ho/ At ~ 0.28 (~ 17 T) leading to a resonance in (^^(h) (Fig. E]). 

Note that the modified Va has three correction contributions: the first is due to the renormalised a-o coupling of elastic 
strains and internal displacements (see also Eq. while the second is due the direct magnetoelastic coupling and 

the third one due to indirect coupling to 4f quadrupoles via the optical modes. In the adiabatic limit when Uo <IC ^im 
the third term has the same form as the second and just adds to an effective magnetoelastic coupling = g^ + {Ego)^- 
This was already noticed e.g. in Ref. d where it was remarked that the contribution of the a-o coupling leads to an 
anomalously large g'^ as obtained from the temperature dependence of the sound velocity. It should be noted however 
that the adiabatic condition does not hold for TGG where the optical mode frequency is quite large as compared to 
the important CEF transition energies (Sect.|Vl. In this case the exact formula for Va{T) given in the above equation 
has to be used. 

A remark on the solution of Eq. p3p is appropriate: In the nonresonant case where ^ '^fm renormalised Wo 
and may be obtained by inserting the unperturbed uJo on the right hand side of Eqs. p3|34|) . Due to the change 
of CEF levels with applied fields the resonance condition ui'^ ~ Af^ for some excitation Ei' Em' may be achieved 
at a special field ho. In this case optical phonons and the CEF excitation will form 'mixed modes' which have no 
longer pure phononic or CEF excitation character. Then Eq. (|33|) has to be solved selfconsistently. Neglecting the 
nonresonant contributions of all other excitations and defining A[,^, — Ai'rn'iPm' ~Pi')S'^^" > one obtains: 

= i(c.2 + Af,„,,) + a.m' - Hm'f + ~9WoH' ra'] ^ (35) 

where a;'m' = sign{uj^ — Af,^,). Passing through the anti-crossing point of CEF excitation Aiim'{h) and phonon 
frequency lUq the frequency of the phonon-type branch changes only by a small finite amount given by {5o/ujo)T^n = 
2{SIj^' gl/iJo)^ < 1 as long as gl/ujo < 1- 

V. APPLICATION TO TGG: CEF STATES AND THEIR ZEEMAN SPLITTING 

A prerequisite to calculate the quadrupolar susceptibilities that enter the Faraday rotation angle is the knowledge 
of CEF states and their Zeeman splitting in a magnetic field. The O]^ space group symmetry of TGG leads to a D2 
site symmetry of Tb'^+ (J=6) non-Kramers ions. The corresponding CEF level scheme has been given in Ref. [T^. It 
may be approximated by a cubic doublet-triplet-singlet level scheme r3(Ad = K) - r5(At = 57.3 K) - r2 (A^ = 
190 K) [la] obtained for parameters W=-24.5 and x=0.8525 from the tables of Ref. \VA Higher lying CEF states are 
neglected. 

As far as eigenstates jFJ^} in a magnetic field are concerned it is important to note that the linearly split triplet 
states |f 5' ) = jFg' ) and the field- independent doublet component |f |) = |r§) have the unperturbed wave functions 
(Fig. [1]) whereas the split-off ground state is mixed with a triplet component according to 

\r\^\T\) + '^\Tl) (36) 
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The inclusion of the orthorhombic (D2) CEF term for general field strength requires numerical calculations similar 
to Ref. m. Now we need the matrix elements of the quadrupolar operators O^z and Oyz of Eq. ^ in the basis of 
Zeeman split CEF eigenstates |r"). They only have to be considered in the subspace S3 = {|r§), jEg), jEg)} where 
they have non-zero matrix elements. The latter are given by 

mil = {rl\o,z\rl) = -i{Tl\Oyz\Tl) 

mi3 = {rl\o,z\rl) ^ i{rl\Oyz\Tl) (37) 

Using the (anti-) symmetrized matrices A'^^, S]^ from Eqs. (|A3IA4p in appendix |X] this leads to explicit quadrupolar 
susceptibilities for the TGG CEF states: 

{{d...Az))'^ = {{dyzdyz))'^ = ml^R-{^){p\~pl)+ml^R+{u){pl~pl) 

{{d.Ajz))l = -{{OyAz)): = ml,Ri^{uj)ipl-pl) - m%R+{io){pl-pl) (38) 

Note that for T the CEF occupation differences {pi ~ Pm) approach unity. Furthermore we defined 

We have shown before that because of the different frequency dependences (t>'^ ^ 0^ for acoustic frequencies. Therefore 
(f)' ~ and we will discuss only the latter in the following. For T — >■ we can write, using Eq. (|A6[) 

Fo{uJo,H) = hU{d:,zdyz)y^^ = ~K^tml[mlR-{uj,,H) - mlR+{ujo,H)] (40) 
= ~Kl(l + pao[mlRd{0,H)+mlR+{0,H)]) 



with = ( ) (^5o)^ and p^o — • In TGG we have < at zero field. However increases linearly with 

field strength (wfg < 0) and will eventually cross the optical phonon frequency (Fig. [T]). Then the Faraday rotation 
becomes large due to an optical phonon resonance in R~ {uio, H). Similarly A^ decreases with increasing field and 
vanishes when the lowest triplet component crosses the ground state. In this case the effective a-o coupling strength 
diverges leading to an additional field dependence in the Faraday rotation caused by i?^(0,iJ). To avoid singular 
behaviour in the resonance case a finite linewidth of phonons (F) or CEF excitations (F^ ) should be introduced. This 
can easily be done by replacing i?^ ^(w, H) by the averaged quantities R^ /j(a;, H) which contains the finite linewidths 
(see appendix [B]) . This leads to the final formula for the Faraday rotation in TGG which will be used for numerical 
calculation: 



' (Ar'-C^2)2+2F2(A,-2+c^2)+r4 (A+2 -t^2)2+2r2(A+2+C^2)+r4 ^ ' 



Here the pre- factor kf — k^/ k^ describes the field dependence of the effective a-o coupling strength. It is determined 
by Pao = {'^nT-Q9^ / E) which gives a measure of the magnetoelastic corrections relative to the bare a-o coupling E. We 
obtain from the last of Eqs. (P7|) : 



k^H) = 1 + Pao {ml ^-T_^^-2 + ^+ A+2 + r+2 ) (42) 



If we ignore these corrections then k^{H) = 1 in Eq. (|4T|). Furthermore the scale of the rotation angle is set through 
^ uja — 2faf which is explicitly given by 

2 {rnqkoY 



A \Va 



(43) 



Note that the Faraday rotation is no longer singular for finite line widths F and F^. However it develops a resonant 
behaviour as function of field which is determined by two effects: The optical phonon resonance with the singlet- 
triplet excitation at A^{H) — Uo and the field dependence of the effective a-o coupling constant in Eq. when 
Af{H) = 0. In TGG these conditions are fulfilled for approximately the same field iJ ~ 17 - 20 T (Fig. [T]) . The 
above treatment may easily be generalised to include coupling to more than one optical phonon by summing over 
individual n-th phonon contributions like Eq. (j4ip each characterised by a frequency Wo"^ and an effective a-o coupling 
strength Kq""* as well as pi"''. 
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h/A h/A 

FIG. 2; Faraday rotation angle (jt'o/vo as function of field strength {ujo = 2.5At,r = 0.2At). Left (a): The optical phonon 
resonance appears around = ujo (crossing with lowest optical phonon mode, Fig. [!}. Full curve: ttiq = 0. Dash-dotted 
line: mg = 0.2. The field dependence of quadrupolar matrix elements enhances the Faraday rotation. The bare a-o couphng 
{pao = 0) is used in Eq. (|42|) Right (b): Rotation angle for renormalised a-o coupling (pao = 0.05, dash-dotted line) in comparison 
with bare case {pao ~ 0, full line). The enhancement factor k^{H) (dotted line) leads to an additional field dependence due to 
level crossing A^ = (Fig. [l| which is close to the optical phonon resonance A^" — ujo- Here a CEF line width Ft — 0.1 At is 
used. 



VI. NUMERICAL RESULTS AND DISCUSSION 



Due to the low D2 site symmetry the precise form of CEF states is not known with certainty. The cubic 
approximation proposed in Ref. [l6| using Ref. [l3 is problematic because it leads to a linear Zeeman splitting 
of Fs triplet which is far too small as compared to the calculation using the full D2 CEF potential performed in 
Ref. fis'l. For a calculation of the Faraday rotation angle according to Eq. (|32l) within the latter approach it would 
be necessary to obtain the field dependent wave functions and energy levels numerically. To gain a principal insight 
a more simplified empirical approach was used here which allows to obtain closed expressions for the rotation angle 
in Eq. (|4ip . We kept the simplicity of the cubic level scheme and wave functions following Ref. and treated the 
dipolar and quadrupolar matrix elements as empirical parameters. The former may be fixed to obtain the qualitative 
level splitting and the latter enter only as scale factors that determine the strength of the resonance behaviour of (l)'{h). 

The simplified CEF level scheme is shown in Fig [TJ The essential states for the Faraday rotation are the split-off 
r| ground state component and the two linearly split Fs-states. Quadrupolar transitions between them lead to the 
rotation angle according to Eq. (HIT) . While the transition energy decreases, (wgg < 0) increases with field 
strength. The field ho where the latter crosses the optical mode at ojo — 2.5At is indicated by an arrow. 

The Faraday rotation angle is shown in Fig. [2l It vanishes for /i — which is in contrast to the Faraday rotation 
for ferromagnets where the time reversal symmetry is already broken by the spontaneaous magnetisation. Under 
such condition the rotation angle is non-zero even at zero external field 0, At the field where A^{ho) — oJo the 
rotation angle in Fig. [2^ becomes resonant and changes sign when the field sweeps across. This is enhanced by the 
ground state level crossing when A'^ = at a slightly larger field which influences the effective renormalised a-o 
coupling (Fig. [^b)- 

The qualitative behaviour of Faraday rotation in Fig. [2] is very similar to the one observed in ferromagnets like 
YsFesOis (YIG) 0. The shape of the resonance depends on the phonon and CEF excitation linewidths and its 
amplitude is enhanced by the opposite field dependence of the quadrupolar matrix elements m± (Eq. (|A6p) as well 
as a-o coupling renormalisation pao (Fig.[2]). 

In an experiment the rotation angle (j)'{h) cannot be observed directly but is infcrcd from the intensity modulation 

I{h;zo) = A^cos^[cl)'{h)zo] (44) 

of the x-polarisation component of the propagating sound wave as it is detected at the end of the sample of length 
Zq- The typical intensity oscillations are shown in Fig. |31 When the field sweeps through the resonance region the 
Faraday rotation becomes large and therefore the oscillations become rapid. This corresponds qualitatively to the 
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h/A, 

FIG. 3: Intensity of the soundwave component (Eq. (|44|l 'l || x. The rapid oscillations occur for fields with large Faraday rotation 
angle (p'o (Fig. [2]). Parameters are F — 0.2At and mg = pat ~ and (jt'o ~ zo = 1 {zo ~ sample length). 

experimental observations where the resonance is around _ff = 17 — 20 T. In the latter the amplitude in the resonance 
region will also be strongly damped which is not described by the present treatment. Therefore instead of measuring 
the amplitude or intensity one may also measure the complementary damping which is large in the resonance region 
and also exhibits the oscillations due to Faraday rotation. If one has more optical phonons with different Wo"^ which 
couple to the acoustic modes the region of rapid Faraday rotation may exist in a more extended field region. The 
damping of the Faraday amplitude in the optical phonon resonance region can in principle be obtained by using the 
propagators with complex optical phonon poles at ujo + *r and solving for complex wave numbers QLtQr- 

Finally we mention that the present theory also explains the Faraday rotation in the case of paramagnetic CeAl2 
where the effect and its linear frequency dependence has first been observed. In this compound a strong a-o coupling 
may be concluded from the large temperature effects which point to an enhanced effective magnetoelastic coupling 
due to optical phonons 0. Indeed in CeAl2 the Ce atoms form a diamond type sublattice which has acoustic (Ti^) 
and lowest optical (T2g) phonons with the same E (C4i,) symmetry for wave vector q=(q,0,0). This enables the a-o 
coupling as in TGG. In fact inelastic neutron scattering experiments |19| found that a-o coupling in this compound 
is strong enough to cause a large anticrossing effect of acoustic and optical phonons for q about halfway to the zone 
boundary. Therefore the Faraday rotation model discussed here is also perfectly applicable for this compound. It 
should be mentioned however that the perturbative approach which was used here to solve Dyson's equations in an 
external field may not be adequate for CeAl2 due to bound state formation of optical phonons and CEF excitations 
caused by large magnetoelastic coupling. [20j 



VII. CONCLUSION AND OUTLOOK 



The origin of acoustic Faraday rotation in paramagnetic 4f compounds has remained mysterious for a considerable 
time. Simple magnetoelastic theories consistently predicted that the rotation angle per unit length should behave like 
(f>' ^ whereas experimental results ^ [l^ invariably have shown that (f>' ^ uja in the case where different sound 
frequencies uja have actually been used. 

In this work we have found the origin of this discrepancy and proposed a theoretical model for TGG which 
resolves the issue. In non-Bravais lattices acoustic and optical phonons may have the same symmetry with respect 
of the group of the wave vector q. In this case long wave length acoustic and internal optical displacements are 
intrinsically coupled both via background elasticity and the magnetoelastic coupling to CEF states. In such a case 
the acoustic Faraday rotation cannot be separated from the field splitting of optical phonons in a magnetic field. 
This leads to indirect contributions (j)'^ to the former which are proportional to the effective a-o coupling squared. 
These contributions were shown to lead to a behaviour 0^ ~ iOa in agreement with experimental observation. They 
dominate the purely acoustic contribution (j)'^ ^ uf^ due to the smallness of sound frequencies. Furthermore it was 
found that in TGG the crossing of CEF triplet levels with the lowest optical phonons and with the lowest doublet 
ground state component leads to resonant behaviour in the Faraday rotation angle for a field strength H ^17 — 20 
T. This has indeed been found in static and pulsed field experiments p^ . 
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We believe the mechanism proposed here may be of more general validity. Certainly it is present in CeAl2 where 
the paramagnetic acoustic Faraday rotation was first found. Indeed the diamond Ce sublattice of this compound 
ensures that acoustic and low-lying optical phonons have the same symmetry and are therefore strongly coupled for 
propagation along the cubic axis. Because Ce^+ is a Kramers ion as opposed to Tb'^^ its CEF level scheme and 
therefore details of the theory are, however, quite different. 

Finally we would like to propose that the field dependence of optical phonons in TGG which approximately extend 
from 90 — 900cm~^ should be investigated by IR absorption and Raman scattering in high fields. According to our 
analysis there should be a field induced splitting of some of the doubly degenerate optical phonon modes which is 
indirectly responsible for the acoustic Faraday rotation. The observation of this splitting, in particular for the low 
frequency optical modes would be a direct support for the theory presented here. 
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Appendix A 

Here we derive the CEF energies and quadrupolar matrix elements of the field-split TGG (cubic) level scheme. 
Using the explicit form of CEF states given in the Lea, Leask and Wolf (LLW) tables [13] , the real matrix elements 
of the dipolar operator in the Zeeman term are given by 



55 



= (ri|j,|ri) = ~(ri|j,|ri) 



m 



Is = {Tl\J.\Tl) (Al) 



to: 



25 



= (r^k.|r2) 



All other matrix elements vanish. The dipolar matrix elements calculated for the cubic CEF states are unreliable 
due to the admixture of other states by the D2 part. Therefore we treat them as empirical parameters to reproduce 
the qualitative Zeeman splitting of triplet states. For magnetic field applied along a cubic axis the Zeeman energy is 
Hz = —hJz- Up to second order in h the CEF level splitting and shifts of doublet, triplet and singlet are obtained as 



At 



El = Ad = 



^3 

pi, 3 



E'^^' = AtTml^h (A2) 



Ei = At = At 



E2 = A, = A, 



A, - At 



This leads to the symmetrized and antisymmetrized quadrupolar matrices of Eq. ([3T|) in the 5*3 — {[f |), iFj), \Tl)} 
subspace according to 

Sxz,xz — Syz^yz = I nT-ll I ; Axz.xz — -Ayz.yz — (■'^3) 





^xz,yz — ~-^yz,xz — I ^^11 | J Sxz,yz — Syz^xz — (-^4) 



TO? 1 TO?3 \ 




; 


A —A — n 

-^xz.xz — -^yz.yz — ^ 


/ 




-n^ii \ 







1 Sxz,yz — Syz^xz — 


/ 





The excitation energies associated with non-zero matrix elements between ground state and two excited triplet com- 
ponents are given by (Fig. [l}: 

At" = El-El=At[l-5 + e^] 

A+ = El-El^At[l + 5 + e^] (A5) 
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where S — and e — "^^^ . Furthermore the field dependent quadrupolar matrix elements may be written as 

rriQ 

mil = m-Q(l -<^) = mnrh- 

TOi3 = -mQ(H e) = -mQm+ (A6) 

where mg is the zero-field singlet-triplet quadrupolar matrix element which may be expressed by the coefficients of 
the singlet-triplet wave functions. It enters into the proper scale Eq. ((43|) of the rotation angle. 



Appendix B 

Finally we give the frequency dependences of diagonal and off-diagonal quadrupolar susceptibilities which have been 
averaged over the phonon and CEF spectral line shapes. For the diagonal susceptibilities we introduce a linewidth 
for the relevant triplet excitations . Then we have {uja 0) 



oo 



Tf/n 2 A 2Af 



R'-^i-a^H) - J_JA^^-^^^^-^^—, - (^p:^ (Bl) 

Likewise for the non-diagonal susceptibilities we introduce a linewidth F for the optical phonon with frequency Wo 
leading to {a — ±): 



R - / du: _ ^^^2 ^ p2 ^a2 _ ^2 ~ (A«2 _ ^2)2 + 2T^{/\f + Uol) + F* ^^^^ 



These expressions are used in Eq. (|40l) to obtain the result in Eqs. (|41l42p 
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